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Abstract 

Recent results show that when non-linear electrodynamics is considered the 
no-scalar-hair theorems in the scalar-tensor theories (STT) of gravity, which are 
valid for the cases of neutral black holes and charged black holes in the Maxwell 
electrodynamics, can be circumvented [U [2] . What is even more, in the present 
work, we find new non-unique, numerical solutions describing charged black holes 
coupled to non-linear electrodynamics in a special class of scalar-tensor theories. 
One of the phases has a trivial scalar field and coincides with the corresponding 
solution in General Relativity. The other phases that we find are characterized by 
the value of the scalar field charge. The causal structure and some aspects of the 
stability of the solutions have also been studied. For the scalar-tensor theories 
considered, the black holes have a single, non-degenerate horizon, i.e., their causal 
structure resembles that of the Schwarzschild black hole. The thermodynamic 
analysis of the stability of the solutions indicates that a phase transition may 
occur. 



1 Introduction 

Among the most natural generalizations of General Relativity (GR) are the scalar- 
tensor theories of gravity in which a single or multiple fundamental scalar fields are 
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present as a possible remnant of a fundamental unified theory like string theory or 
higher dimensional gravity theories [3] . Different modifications of scalar-tensor theories 
are attracting much interest also in cosmology and astrophysics. 

Possible deviations from GR, especially the existence of new effects, due to the 
presence of the scalar field would be of considerable interest. According to the no- 
scalar-hair conjecture, the black-hole solutions in the STT coincide with the solutions 
from GR. No-scalar-hair theorems treating the cases of static, spherically symmetric, 
asymptotically flat, electrically neutral black holes and charged black holes in the 
Maxwell electrodynamics have been proved for a large class of scalar-tensor theories 
[H EJ [6]. The scalar field in these cases is constant, and thus trivial, if one demands 
that the essential singularity at the center of symmetry is hidden in a regular event 
horizon. 

There are no no- scalar-hair theorems, however, in the non-linear electrodynamics 
(NLED). It was first introduced by Born and Infeld in 1934 to obtain finite energy 
density model for the electron [7]. For more information on recent interests in non- 
linear lagrangians of electrodynamics please see [8], also [9] [30] and references cited 
therein for gravitational aspects of NLED. 

In the case of NLED, the energy-momentum tensor of the electromagnetic field 
has a non-vanishing trace which is non-trivially coupled to the scalar field. Hence, 
the electro-magnetic field acts as a source of the scalar field and allows the existence 
of asymptotically flat, hairy black holes. Such solutions have been found recently in 
different non-linear electrodynamics [H [2] . These solutions are hairy in a sense that 
the scalar field is not trivial. That hair, however, is secondary since the solutions are 
determined uniquely by the values of their magnetic charge and mass, and the value of 
the scalar field at infinity. For other solutions describing asymptotically flat black holes 
with scalar hair (in which, however, the potential of the scalar field is not positively 
semi-definite) we refer the reader to |31j . 

Another interesting effect that occurs when NLED is considered is the presence of 
multiple black-hole phases in a certain class of STT. Presence of non-unique solutions 
in STT has been found by Damour and Esposito-Farese in their works on neutron 
stars in the STT. One of the most interesting effects that occur in these solutions is 
the spontaneous scalarization, a scalar field analogue of the spontaneous magnetization 
of ferromagnets [321 [331 [31]. In the present work we consider the same STT and find 
multiple-phase black- hole solutions numerically. Unlike the situation in [UI2], they are 
not determined in a unique way by the values of their magnetic charge and mass and 
the value of the scalar field at infinity. The phase diagram of the present solutions is 
reminiscent of the phase transition between caged black holes and black strings in higher 
dimensions [35l 136] . The preliminary analysis gives us a reason to suppose that in the 
system we study a phase transition might occur. The analysis of the thermodynamics 
of the solutions, however, is not sufficient to make a conclusion about the stability of 
the solutions and about the existence of phase transitions. The full examination of the 
problem requires a perturbative analysis. 
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2 Basic equations and qualitative investigation 

The theory we consider is presented [H [2] both in the Jordan and in the Einstein 
frames. In the Einstein frame the action of the theory we consider is 

S = -tV / d 4 x^-g [U - 2g^d fi( pd u i P - 4V(cp) + AA\ V )L(X, Y 

V(ip) is the potential of the scalar field, L(X) is the Lagrangian of the electromagnetic 
field, A 4 (tp) is the function which determines the coupling between the electromagnetic 
field and the scalar field. We also have that 

X = ^ELf^^f^ Y = (*Fy» (1) 

where "★" stands for the Hodge dual with respect to the Einstein frame metric g^ v . 
The type of the STT is determined by the explicit choice of the functions V(ip) and 
A 4 (f). The field equations obtained from action (CQ) take the following form 

= 2d^d v p + 2V(ip)g ia , - 2d x L(X, Y) (f^F? - ^g^F^F ^ 
-2A 4 &) [L(X, Y) - Yd Y L(X, Y)] 9fMV , 

V, [d x L(X, Y)F^ + d Y L(X, Y)(*Fy»] = 0, (2) 

V M W = ^r^- ~ M<f)A\<P) [L(X, Y) - Xd x L(X, Y) - Yd Y L(X, Y)} , 

dip 

where a((p) = dln ^ y ' > . 

In what follows the truncated^ Born-Infeld electrodynamics described by the La- 
grangian 



L BI (X) = 2b\l- ) Jl + jj (3) 
will be considered. And V(<p) will be equal to zero. 

The anzats for metric of a static, spherically symmetric space-time can be taken in 
the following form 

ds 2 = g^dx^dx" = -f(r)e~ 2S{r) dt 2 + — — + r 2 (d6 2 + sin 2 6d(f) 2 ) . (4) 

/(?") 

Since the Born-Infeld NLED is invariant under electric- magnetic duality rotations we 
will study only the magnetically charged black holes for which the electromagnetic field 
is given by 

F = P sin 9d6 A d<f> (5) 
and the magnetic charge is denoted by P. 



^^Here we consider the pure magnetic case for which Y = 0. 
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The field equations reduce to the following coupled system of ordinary differential 
equations: 



dS ( dip\ 

dr \dr J ' 



dm 
dr 



H%) - Al[mx) 



— ( r 2 f— \ 
dr \ dr J 



where 



r 2 <^ -Aa{p)A\p) [L(X) - Xd x L{X)} - rf 



2m 




(6) 
(7) 

(8) 



/ = 1 



and X reduces to: 



X 



A~\p) P 2 



(9) 

We will be searching for solutions which have a regular horizon on which the scalar 
field p does not diverge. The regularity of the transition between the Einstein and the 
Jordan conformal frames requires the following restrictions on the coupling function 
< A(ip) < oo for r > th, where r# is the radius of the horizon, and we will consider 
STT for which it is satisfied. The diversity and the properties of the solutions depend 
strongly on the choice of the functions A(<p) (respectively on «(</?)). Even when the 
functions A(p) are chosen to satisfy the experimental constraints (see, for example, 
|33j ) we are left with infinitely many possibilities. That is why some restrictions on the 
functions A((p) representing the STT should be imposed . In the present work we will 
consider only theories for which ip a(ip) is non-negative or non-positive for all values of 
(p. As we prove below, theories for which ipa(ip) > for all values of cp do not admit 
(asymptotically flat) black hole solutions with non-trivial scalar field. In theories with 
ipa(ip) < black holes (if they exist) have a single non-degenerate event horizon when 
a ((f) = admits only the solution <p = 0. 
Using the following equation 



d 
dr 



dr I 



Ar 2 e- s a(ip)A\ip) [Xd x L(X) - L(X)} 



(10) 



which is another form of equation 
© holds 

Xd x L{X) 



and the fact that for the Born-Infeld Lagrangian 



L(X) > 



we can draw some conclusions about the general properties of the solutions. 

Let us multiply equation ( fTUl) by ip and then integrate it in the interval r e [r#, oo) 
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where we denote the radius of the outer horizon (the event horizon) with r# 



d 
dr- 



rii 



4 / r 2 e~ s ipa(tp)A\ip) [Xd x L(X) - L(X)) dr, 



(12) 



in 



and after integrating by parts we get 



,• I -8 2 f d V\ ( -8 2, dip 

urn I ipe r j— J — I ipe r j — 



dr 



dr 



lim T>ip — 



r=rjj 



in 



-5 2 f f d( P 



dr 



rH 



-8 2 f ( d( P 

\ dr 



dr 



rV V a((^)^ 4 (y) - L(X)] rfr. (13) 



'7/ 



In the second line we take advantage of the fact that we are looking for asymptotically 
flat, black hole solutions so }'{r H ) = 0, lim^oo f(r) = 1, lim^oo 8(r) = and for the 
asymptotic value of the scalar field we impose lim^oo tp(r) = ip^, where tp^ = 0. The 
constant V denotes the scalar charge which is defined as 



D = — lim r 2 -^. 

r^oo dr 



(14) 



Since ip is vanishing at infinity we finally get 



oo 




/ 


e- 5 r 2 f 


r H 






oo 




= ifr 







\dr ) 



dr 



4 / r 2 e- 5 <pa(ip)A 4 (<p) [Xd x L(X) - L{X)\ dr. 



(15) 



The total sign of the left-hand side of (fl5l) is negative for black hole solutions with 
nontrivial ip. The sign of the integral on the right-hand side depends on the sign of 
(pa(ip). If <pa.(<p>) > a contradiction is reached so the assumption for the existence 
of asymptotically flat black holes with non-trivial scalar field in this theory is wrong. 
Asymptotically flat black-hole solutions with nontrivial scalar field exist in scalar-tensor 
theories for which (pa((p) < 0. Through a similar examination we prove that these 
black holes have a single, non-degenerate horizon. Let us admit that more than one 
horizon exists. Then we multiply equation fflUj) by <p again and integrate it by parts 
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in the interval r G [r_,r + ] where we denote the first inner horizon and the outermost, 
non-degenerate horizon with r_ and r + , respectively 




\ dr 




T— 



C rf \Tr 



2 



dr 



4 / rVV(#V) [Xd x L{X) - L(X)} dr < 0. 



(16) 



r~ 



Having in mind that /(r_) = = f(r+) and that / < in the interval r G [r_, r+] we 
reach a contradiction, which means that the admission is incorrect. 

Now, we only have to prove the non-existence of extremal black holes (black holes 
with degenerate event horizon). For the scalar-tensor theories we consider a((p) turns 
to zero only when tp = 0. Let us admit that an extremal black hole with non-trivial 
scalar field exists. In this case, the left-hand side of (ITU1) is equal to zero. The right- 
hand side is equal to zero only when = 0, where tpjj is the value of the scalar 
field on the horizon, which means that tp H = 0. We also require that tpoo = 0, where 
ipoo is the value of the scalar field at the spacial infinity. In this situation, the only 
possibility to have a solution with non-trivial scalar field is the scalar field to have at 
least one extremum. So let us integrate equation (imp in the interval r G [r# , r e ] , where 
r e is the point of the leftmost (the one which is nearest to the event horizon) extremum 
of tp which is on the right of the event horizon 



Since tp ^ in the interval (r/r,r e ], the sign of a(<p) also does not change in this 
interval. Hence, the integral on the right-hand side of (ITT)) is non-zero and has a fixed 
sign which depends on the sign of cx(<p)- The contradiction we reach means that no 
extremal solutions with non-trivial scalar field can exist. 

To sum up, we can say that if a black hole exists it will have a single horizon, i.e., 
its causal structure will be Schwarzschild-lik^]. In both conformal frames, inside the 
event horizon a space-like singularity is hidden. 

Finishing this section it is worth noting that the differential equations system (JBT 
IHD is invariant under the rigid rescaling r — > Ar, m — > Am, P — > XP and b — > \~ 2 b 
where < A < oo. Therefore, given a solution to ([6H8]) with one set of physical pa- 
rameters (rh, M, P, b, V, T), the rigid rescaling produces new solutions with parameters 
(Arft, AM, XP, X~ 2 b, XV, A _1 T). Here T denotes the temperature of the horizon. 

2 Another class of scalar-tensor theories which admit black holes of the Schwarzschild type are those 
with negative function /3(<p) — da }^ for all values of ip and a((^oo)=0. This can be shown by a method 
similar to that presented above. One can also show that scalar-tensor theories with /3(<p) > for all 
values of ip do not admit asymptotically flat black holes with nontrivial scalar field. 




(17) 
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3 Numerical results 



The nonlinear system (jSJ-dHI) is inextricably coupled and the event horizon r# is a priori 
unknown boundary. In order to be solved numerically, it is recast as a equivalent first 
order system of ordinary differential equations. Following the physical assumptions of 
the matter under consideration the asymptotic boundary conditions are set, i.e., 

lim m(r) = M (M is the mass of the black hole in the Einstein frame) , 

r— >oo 

lim S(r) = lim (p(r) = 0. 

■r— >oo r^oo 

At the horizon both the relationship 

f(r H ) = 

and the regularization condition 

= {Aa^)A\ V )[Xd x L(X) -L(X)]}\ r=rH 

concerning the spectral quantity r H must be held. For the treating the above posed 
boundary- value problem (BVP) the Continuous Analog of Newton Method (see, for 
example |37j,[38j,[20]) is used. After an appropriate linearization the original BVP is 
rendered to solving a vector two-point BVP. On a discrete level sparse (almost diagonal) 
linear algebraic systems with regard to increments of sought functions 5(r), m(r), and 
ip(r) have to be inverted. 

For our numerical solutions we have considered the STT studied by Damour and 
Esposito-Farese in their works on neutron stars in the STT of gravity. In this particular 
theory, the coupling function has the following form 

A(<p) = e> 2 , (18) 

where (3 is a negative constant. Observational data from binary-pulsar and solar- 
system experiments restricts the admissible values of (3. When a(<p = 0) = 0, as it is 
in our case, the coupling constant should be (3 > — 5 (see [33] for more details). We 
have studied the solutions for values of the parameters which are in agreement with 
the current observations but also for such values that are out the admissible interval 
since the later have qualitatively different behavior from the former which makes them 
interesting for the theory. 

For this coupling function the field equations possess the discrete symmetry ip — > 
—(p. Let us also note that every general relativistic solution is a solution to the scalar 
tensor theory under consideration with ip = 0. 

A thorough study of the phase space would be difficult due to the large number of 
parameters (/5, b, P, M) in the problem. So in order to illustrate the general behavior of 
the obtained solutions we give several representative figures considering several values 
of (3, b and P and varying the mass M of the black hole. For the cases presented here 
6 = 0.01; 0.2. 



df dip 
dr dr 
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3.1 General description of the phase space 



Even if the values of all four parameters /?, b, P, M and the boundary conditions are 
fixed the solutions of ([ME]) are not uniquely determined, i.e. the solutions are not 
unique. An additional parameter should be used for labeling of the different solutions. 
One natural choic^f) would be to label the different solution by the value of the scalar 
charge T>, which here unlike the cases in [TJ [2] is independent. 

The global structure of the phase space changes with the variation of (3. For f3 > 
/5crit, where (3 crit ~ —14.9 when b = 0.01 and P = 1.0, the M—V phase diagram consists 
of three branches, while for (5 < j3 CT - lt the number of branches is five. The change of the 
qualitative structure of that phase diagram with the variation of (5 is given in Figure 
(CQ) . The different cases have been studied in more details in the succeeding subsections. 
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Figure 1: The M — V relation for several different values of /3, b — 0.01 and P = 1.0. 
For (3 > Pent the Middle branch disappears. 



3.2 Cases with three solutions, j3 = —4.0 

In the left panel of Figure (T2]) the dependence T>(M) of the scalar charge on the mass 
for P = —4.0 is presented. There are two special points to be considered, namely A 
and E. For P = 1.0 the points A and E lie at masses M A ps 0.794 and M E w 0.46. For 
values of the black hole mass in the interval M G (Mb, Ma) three solutions co-exist. 
We call them Outer and Trivial. The equations posses a discrete symmetry ip — > —ip 
as a result of which the Outer solution has a mirror image with respect to the abscissa, 
whose scalar field charge T> has an opposite sign. The symmetric solutions we will 
denote as Outer ± where the indices + and — refer to the sign of the scalar field charge 
of the solutions. The radius and the temperature of the event horizon, are the same 
for both solutions in the couple so when we comment on them will usually omit the + 
and — indices. 

3 An alternative choice would be the value of the scalar field on the horizon 
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The Trivial brunch has scalar field ip = and represents the Einstein-Born-Infeld 
solution in GR. For masses lower than Me only the solutions which we call Outer exist 
and for M > M4 only the Trivial solution remains. 

The radii of the black holes are shown in the right panel of Figure ([2]). An object 
with zero radius of the event horizon (a naked singularity) is reached for a finite value 
of M for all three solutions - the Outer± and the Trivial. 

The inverse temperature T _1 of the solutions is presented in Figure ([3]). 

3.3 Cases with five solutions, (3 = —50.0 
3.3.1 Solutions with Born-Infeld parameter b=0.01 
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Figure 5: The M — r# relation for the same values of the magnetic charge as in Figure 
03]). A magnification in the encircled region is presented in Figure (jBJ). 

For (3 = —50.0 and b = 0.01 the T>(M) dependence is given in Figure (jlj). Here, the 
special points to be considered are four, A, B± and E. For P — 1.5 the points A and 
B± lie at masses M A w 2.023 and M B± = Mb ~ 2.37 while in the case of P = 3.0 we 
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have Ma ~ 3.60 and M E± = Mb ~ 4.53. As it can be seen, for values of the black hole 
mass in the interval M G (Ma, Mb) five solutions co-exist. They can be separated in 
three groups which we name Outer, Middle and Trivial. The symmetric solutions we 
will denote as Outer ± and Middle± with the same convention for the + and — indices 
as in the previous case. Both points B + and S_ are projected on the same point on 
the M - r H and M - T _1 dia grams, which we denote simply as B. 

The Trivial has scalar field <p = and represents the Einstein-Born-Infeld solution 
in GR. For masses lower than M E ~ 1.00 and M E ~ 2.42 for P = 1.5 and P = 3.0, 
respectively, only the solutions which we call Outer exist and for M > Mb only the 
Trivial solution remains. 

The radii of the black holes are shown in Figures (j3j) and (jSj). Again, for three of 
the solutions - the Outer and the Trivial a naked singularity is reached for a finite 
value of M. The Outer solutions have a larger radius than the Trivial for low masses, 
but as it can be seen in Figure (jSJ), which is a magnification of the encircled region in 
Figure flSD, with the increase of the mass, in point C, the situation changes and the 
black hole with zero scalar charge becomes larger. The approximate position of point 
C is shown also on Figure (j3J) with a dotted vertical line. The Middle solution black 
holes are smaller than the other three for all values of M. 

The inverse temperature T~ x of the solutions is presented in Figure (j7j). With the 
decrease of the mass M the inverse temperature of the Trivial solution passes through 
a local maximum which gets sharper with the increase of the magnetic charge P and 
leads to numerical calculation difficulties. In the limit of vanishing radii of the black 
holes their inverse temperature decreases steeply. 

3.3.2 Solutions with Born-Infeld parameter b=0.2 




0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 



-0,8-1 - " M 

Figure 8: The M — T> and the M — r^ relations, for b = 0.2 and P = 1.0. 

Here we present an example of solutions for values of the parameters for which the 
Trivial solution reaches an extremal black hole instead of a naked singularity with the 
decrease of the mass M. In Figure (ED the dependence T>(M) of the scalar charge on 
the mass and the radii of the black holes are presented. The points A and B± lie at 
masses Ma ~ 1.13 and M E± = Mb ~ 1.54. Again, for values of the black hole mass 
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in the interval M G (M4, Mb) five solutions co-exist. The Trivial reaches an extremal 
black hole at Me ~ 0.93 and this can be seen on the M — r# diagram. 

In Figure a magnification of the encircled region from Figure (jHJ) and the inverse 
temperature are shown. The point C is once again indicated with a vertical line. 
With the approaching of the extremal black hole the inverse temperature T _1 rises 
unboundedly. 



4 Thermodynamics 

Black holes have long been known to be thermodynamical systems [39]. The First Law 
(FL) of black hole thermodynamics in the presence of a scalar field has the following 
form [10] 

6M = T6S+V H SP + V6<p 00 , (19) 

where T, S and P are the temperature, the entropy, and the magnetic charge of the 
black hole, respectively, and <£>oo is the value of the scalar field at spacial infinity. 

The quantity \I> conjugate to the magnetic charge is the potential of the magnetic 
field which is given by the following definition 

H„ = d^. (20) 

On the other hand the magnetic field is defined as 

H, = -^G^C, (21) 

where 

G - = - 2 ^P* (22) 

£ = is the Killing vector generating time translations and "*" is the Hodge star 
operator. 
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Since in our case the asymptotic value of scalar field is fixed, the term which contains 
its variation vanishes and the FL reduces to the form it has in GR. 

In the situation when the solutions are not unique a natural question is which of 
them are stable. Certainly, the answer of this question requires linear perturbative 
analysis of our system of coupled differential equations and solving the corresponding 
eigenvalue problem. In certain cases, however, some information on the stability can be 
inferred by using only the equilibrium thermodynamical characteristics of the solutions 
via the so-called "Poincare" or "turning point" method. For a nice discussion of the 
method we refer the reader to [H] and references therein. The advantage of this 
method is its remarkable simplicity. The "turning point" method, however, may hide 
many uncertainties and should be applied with caution. The method consists in the 
following. Consider a system with thermodynamical parameters /i\ The equilibrium 
states (stable or not) are extrema of an appropriate Massieu function S. At equilibrium 
state we can define the conjugate variables Pi(id) = ^f- According to the method the 
change of stability can only occur at turning point£| or bifurcations in the equilibrium 
sequence on the conjugate diagram In the absence of bifurcations the stability 

character changes only when the equilibrium curve meets a turning point and if one 
single point of an equilibrium sequence is shown to be fully stable, then all equilibria 
in the sequence are fully stable up to the first turning point. At the turning point the 
branch with negative slope is always unstable while the branch with a positive slope is 
more stable than the one with negative slope. 

Concerning the application of the "turning point" method to our case we shall 
consider scalar-tensor black holes in the micro-canonical ensemble. In this case the 
corresponding Massieu function is the entropy We will comment on the case 

with five branches first. The conjugate variable is J3m = T' 1 . For (3 < /? crit the 
conjugate diagram M — T _1 is shown on Figures ([7j) and (Q. The uncertainties in our 
case come from that fact that we do not know the full diagram, i.e., whether there are 
other branches and bifurcation points different from point A. Assuming however that 
there are no other bifurcation points and taking into account that there is only one 
turning point B we may conclude that the Middle branch is unstable while the Outer 
branch is more stable. The cusp which appears in point B on the M — Th diagram 
in Figures (jSJ) and (jUJ) also indicates a change in the stability of the solutions. Since 
the Outer branch is the unique solution (up to the discrete degeneracy Outer -±) for 
sufficiently small masses one might accept that the Outer branch is probably stable 
there. Then according to the "turning point" method the Outer branch should be fully 
stable up to the turning point B. The Trivial branch is unique solution for sufficiently 
large masses and as a GR solution is known to be stable^ there [321 [32]. So, assuming 
that Trivial branch is stable for large masses, the "turning point" method asserts that 
the Trivial branch should be stable up to the bifurcation point A. Since the entropy of 
a black hole is proportional to the area of its event horizon, among the three black-hole 

4 The turning points are points where two equilibrium branches merge with a vertical tangent. A 
bifurcation point is point where branching of equilibrium sequences occurs. 
5 We keep the magnetic charge fixed. 

6 In general, the stability of the trivial solution within the framework of GR does not guarantee its 
stability within the "larger" scalar-tensor theory. 
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phases we consider, the one with the biggest radius would have the maximal entropy 
(see Figures (JSD , dSD , (JHJ) and So for M < M c the Outer solutions would be 

thermodynamically favorable and for M > Mq - the Trivial. The Middle solutions 
are thermodynamically unstable for all values of the mass since their radius is smaller 
than the radii of the other three solutions. The point C is a candidate for a point of 
a first order phase transition. Let us stress again that above analysis based on the 
"turning point" method is only suggestive and cannot serve for a definitive solution 
of the stability of the solutions. Reliable analysis of the risen questions will be given 
elsewhere together with the solution of corresponding eigenvalue problem. 

For j3 > (3 CT it the M — T^ 1 diagram is given in Figure ([3]). No turning points can be 
seen on it. As in the previously discussed case, since Outer ± are the only solutions for 
sufficiently small masses one might accept that they are probably stable there. Then 
according to the "turning point" method the Outer branch should be fully stable up to 
the bifurcation point A. Again, assuming that Trivial branch is stable for large enough 
masses one can expect that to the left it should be stable at least up to the bifurcation 
point A. From the right panel of Figure (jSJ) we can see that the Outer solutions have 
larger radius of the event horizon so they would be thermodynamically favorable. 

The thermodynamical stability considerations were made in the Einstein frame. In 
order to transfer the conclusions to the physical, Jordan frame properly, we need to 
clarify the connection between the thermodynamic properties in the two conformal 
frames. The temperature of the event horizon is invariant under conformal transfor- 
mations of the metric that are unity at infinity jUj. The properly defined entropy is 
also the same in both conformal frames. It has been proved that in the Jordan frame 
the entropy of the black hole is not simply one fourth of the horizon area 06] as 
in the Einstein frame and needs to be generalized. The entropy in the Jordan frame is 
defined as 

Sj = ^J<Pxy/=&>jjF($). (23) 
Passing to the Einstein frame we get 

Sj = 4^I ^ X V°^ = S E = S. (24) 

In the last two equations ( 2 'g and ( 2 'g are the determinants of the induced metrics on 
the horizon in the Jordan and in the Einstein frame, respectively. 

The term in the FL (fl9|) connected with the magnetic charge is also preserved under 
the conformal transformations. 

In order for the FL of thermodynamics to be satisfied in the Jordan frame the mass 
should be properly chosen since the Arnowitt-Deser-Misner (ADM) masses in both 
frames are not equivalent. It can be easily shown that in the STT considered the ADM 
mass in the Jordan frame Mj is equal ADM mass in the Einstein frame M 

Mj = M. (25) 

For the Jordan frame, the proper mass in the FL of thermodynamics is the ADM mass 
in the Einstein frame M. Similarly, for boson and fermion stars the proper measure for 
the energy of the system is again the ADM mass in the Einstein frame M. For more 
details on the subject we would refer the reader to the works (13, HHJ H9l [50] . 
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5 Conclusion 



In the present work new numerical solutions describing multiple-phase, charged black 
holes coupled to non-linear electrodynamics in the scalar-tensor theories with massless 
scalar field were found. Since an electric-magnetic duality is present, here only the 
purely magnetically charged case was studied. For the Lagrangian of the non-linear 
electrodynamics the truncated Born-Infeld Lagrangian was chosen and a special STT 
was considered. As a result of the numerical and analytical investigations, some general 
properties of the solutions were found. In the class of STT for which (pa((p) < for all 
values of <p, black holes with a single horizon exist, i.e., their causal structure resembles 
that of the Schwarzschild black hole and is simpler than the corresponding solution in 
GR. In dependance of the black hole mass, different black hole phases, labeled by the 
scalar charge coexist. For sufficiently small masses there are two phases and only one 
phase for sufficiently large masses. In the intermediate range of black hole mass, there 
are three or five coexisting phases. 

The thermo dynamical stability of the solutions was also discussed but the present 
data is not sufficient for a final conclusion. The full understanding of the stability of the 
phases and the dynamics of possible phase transitions requires perturbative analysis. 

Presence of multiple phases is also expected in cases of other sources of gravity 
whose tensor of energy-momentum has a non-vanishing trace, such as: linear and non- 
linear electrodynamics in higher dimensions (or even in 3D gravity), gauge theories 
and so on. 
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